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Description 



Minimizing Computational Complexity 
in Cell-level Noise Characterization 

Background of Invention 

[0001] piQid of the Invention 

[0002] The present invention relates to computer aided design 
(CAD) methodologies used to implement integrated cir- 
cuits, and more specifically to a method and apparatus to 
reduce computational complexity in characterization of 
various parameters representing cross-talk noise of cells 
in a library used in designing integrated circuits. 

[0003] Related Art 

[0004] Libraries are often used in designing integrated circuits 
(ICs). In general, a library contains multiple cells (e.g., 
logic gates, multiplexers, etc.), and a designer selects and 
includes a desired cell in a circuit design. Several libraries 
are often provided, each being designed for a corre- 
sponding combination of manufacturing process, temper- 



ature, voltage (PTV), etc. 

[0005] Page of Cells are often characterized prior to being used 
in the design on ICs. Characterization generally refers to 
the determination of various electrical parameters of in- 
terest for each cell. Such pre-determination enables a 
system (typically, using which an IC is being designed) to 
quickly analyze a proposed design for conformance with 
various criteria (e.g., design rules, timing constraints) as a 
designer designs a desired circuit. 

[0006] Cross-talk noise among nets connecting cells is often of 
concern in the design of ICs, particularly in deep submi- 
cron technologies. Cross-talk noise generally refers to in- 
terference caused by some nets (aggressor) on other nets 
(victims). Of particular concern in IC designs is crosstalk 
noise induced due to parasitic coupling. Such induced 
noise or glitches could cause functional failure by switch- 
ing logic.states. In the case of designs employing dynamic 
logic design styles, this problem is more severe due to in- 
creased noise sensitivities of pre_charged nodes. 

[0007] By analyzing a (proposed) circuit for the effects of such 
cross-talk noise, one may ensure that a later fabricated 
circuit would operate in a desired manner even in the 
presence of cross-talk noise. At least for such a reason, it 



is generally desirable that cells in a library be character- 
ized to determine various related parameters ("noise pa- 
rameters"). 

[0008] One noise parameter of interest is 'noise immunity". Noise 
immunity parameters generally indicate whether a glitch 
of a given width and/or height on an input pin of a cell 
would cause the signal at an output pin (and thus the 
connected nets) to be altered by more than a pre- 
specified threshold (and/or a change of state in the case 
of sequential elements). Often a noise immunity curve 
(NIC) is generated which indicates the minimum height of 
an input glitch that would cause an output glitch exceed- 
ing a pre-specified threshold level/height for a given 
width of an input glitch. 

[0009] Another noise parameter of interest is 'noise propagation' 
(NP), which specifies the characteristics of an output glitch 
in response to a input glitch of a given height and width. 
NP curves (commonly referred to as NP characteristics) are 
often drawn mapping the height (peak) of the input glitch 
to the height of the output glitch. The NP curves would in- 
dicate the effect of an input glitch (on one cell) on the 
subsequent cells connected to the output net of the cell. 

[0010] The characterized NIC and NP curves are used to quickly 



determine failures in response to input glitclies in inte- 
grated circuits. For example, in static noise analysis (SNA), 
the NP characteristics/curves and NICs are used to deter- 
mine the effect of aggressor nets on victim nets. Often the 
analysis is performed without regard to the states of the 
input pins generally because the probability that the 
states of all pins will be known is low. The results of the 
analysis are used as a basis to determine which nets need 
to be redesigned. By pre-characterizing the curves, such 
analysis may be performed quickly. 
[0011] The need for such characterization is increasing with im- 
provements (e.g., from 130 nanometer process technol- 
ogy to 90 nanometer process technology) in fabrication 
technologies, which are enabling packing of more transis- 
tors into the same area. It has been observed that the 
coupling capacitance contributes to a greater degree (as a 
fraction of the total capacitance on a net) with such im- 
provements. Accordingly, an increasing need is believed 
to be presented to characterize the related noise parame- 
ters. 

[0012] One general requirement in characterization of cells for 
the noise parameters often is to reduce the computational 
requirement. Often due to a large number of cells in each 



library and/or a large number of libraries (e.g., to corre- 
spond to different PVT combinations), it may take several 
days or weeks to characterize the noise parameters of the 
cells in cell libraries. By decreasing such total characteri- 
zation time, the total time to design an IC as well as cost 
for characterization, can be reduced. 
[0013] What is therefore needed is a method and apparatus to 
reduce computational complexity in characterization of 
various parameters representing the effect of cross-talk 
noise on cells in a library used in designing integrated cir- 
cuits. 

Brief Description of Drawings 

[0014] jhe present invention will be described with reference to 
the following accompanying drawings. 

[0015] Figure (Fig.)lA is a block diagram of an example circuit 
used to illustrate several aspects of the present invention. 

[0016] Figure IB is a diagram of a cell illustrating various input 
pins and output pins. 

[0017] Figure 2A is a timing diagram of an example glitch repre- 
senting cross-talk noise. 

[0018] Figure 2B is a timing diagram illustrating noise immunity 
curves of a full adder. 

[0019] Figure 2C is a timing diagram illustrating the failure area 



and success area in relation to a noise immunity curve. 

[0020] Figure 3 is a blocl< diagram illustrating various inputs and 
outputs of a full adder. 

[0021] Figure 4 is a flow chart illustrating the manner in which 

computations may be reduced while performing noise im- 
munity analysis according to an aspect of the present in- 
vention. 

[0022] Figure 5 depicts curves illustrating the selection of worst 
case curve when there is no cross-over of the curves fitted 
by a few immunity transition points in an embodiment of 
the present invention. 

[0023] Figure 6 depicts curves fitted from a few immunity transi- 
tion points illustrating the selection of worst case curve 
when top curves are crossing in an embodiment of the 
present invention. 

[0024] Figure 7 depicts curves fitted from a few immunity transi- 
tion points illustrating the presence of more than one 
worst case curve. 

[0025] Figure 8A is a graph illustrating the input glitches applied 
to a buffer. 

[0026] Figure 8B depicts the signals at an output pin of the buffer 

in response to the input glitches of Figure 8A. 
[0027] Figure 9 is a graph illustrating noise propagation (NP) 



curve with the peak values of input glitches shown on X- 
axis and the corresponding peaks of the output responses 
shown on Y-axis. 
[0028] Figure 10 is a flow chart illustrating the manner in which 
computations may be reduced while performing noise 
propagation analysis according to an aspect of the present 
invention. 

[0029] Figure 11 is a flow chart illustrating the manner in which 
immunity transition points on a noise immunity curve may 
be determined using curve fitting technique according to 
an aspect of the present invention. 

[0030] Figure 12 depicts a curve fitted from a few immunity tran- 
sition points illustrating the selection of search range to 
determine immunity transition point according to an as- 
pect of the present invention. 

[0031] Figure 13A is a timing diagram containing input glitches 
and corresponding output glitches while determining NP 
data. 

[0032] Figure 13B is a flowchart illustrating the manner in which 
the search range (to determine an immunity transition 
point) can be determined accurately using NP data ac- 
cording to an aspect of the present invention. 

[0033] Figure 14A is a flowchart illustrating the manner in which 



data related to determined immunity transition points may 
be used to determine the search range accurately while 
generating a noise immunity curve according to an aspect 
of the present invention. 

[0034] Figure 14B depicts curve fitted from a few immunity tran- 
sition points illustrating the selection of a search range. 

[0035] Figure 15 is a flow chart illustrating the manner in which 
various approaches described above can be combined to 
generate a narrow and accurate search range according to 
an aspect of the present invention. 

[0036] Figures 16A and 16B are logical diagrams together illus- 
trating the selection of the search range from the search 
ranges determined in an embodiment of the present in- 
vention. 

[0037] Figure 17 depicts NP curve illustrating the need to place a 
start point in a transition region when determining the 
points of the NP curve. 

[0038] Figure 18 is a flowchart illustrating the manner in which 
the number of simulations can be reduced while generat- 
ing a NP curve according to an aspect of the present in- 
vention. 

[0039] Figure 19 is a flowchart illustrating the manner in which 
the start point may be placed in the transition region for 



the first three NP curves according to an aspect of the 
present invention. 
[0040] Figure 20 depicts first three NP curves corresponding to 
narrowest width, widest width and intermediate (middle) 
width. 

[0041] Figure 21 is a flowchart illustrating the manner in which 
the start point may be placed in the transition region for 
the remaining (after first three) NP curves according to an 
aspect of the present invention. 

[0042] Figure 22 is a flowchart illustrating the manner in which 
the NP curves and the NICs may be determined in an in- 
terleaved manner to reduce computational requirements 
according to an aspect of the present invention. 

[0043] Figure 23 is a timing diagram illustrating the general con- 
sideration in noise immunity characterization for sequen- 
tial elements. 

[0044] Figure 24 is a flowchart illustrating the manner in which 
sweeps can be avoided for some heights of input glitches 
according to an aspect of the present invention. 

[0045] Figure 25 is a timing diagram depicting input glitch and a 
clock signal with large setup and hold times in an embod- 
iment of the present invention. 

[0046] Figure 26 is a flowchart illustrating the manner in which 



noise immunity curves can be generated for sequential el- 
ements according to an aspect of the present invention. 
[0047] In the drawings, like reference numbers generally indicate 
identical, functionally similar, and/or structurally similar 
elements. The drawing in which an element first appears 
is indicated by the leftmost digit(s) in the corresponding 
reference number. 
Detailed Description 

[0048] L Overview 

[0049] An aspect of the present invention enables reduction of 
computational resources in characterizing parameters 
representing cross_talk noise of cells in a library by first 
determining a worst case vector for an input pin and out- 
put pin combination (commonly referred to as an arc), 
which would propagate most noise from the input pin to 
the output pin among all the vectors of interest. A vector 
of interest for an input pin and output pin combination 
represents an input bits combination for other input pins, 
which would cause a bit value transition on the output pin 
if the input bit value is changed on the input pin. 

[0050] An integrated circuit designed using the cell library is then 
analyzed using the parameters data related to the worst 



case vectors. Thus, when a cell is analyzed while receiving 
other (non-worst case) vectors, the parameters data may 
indicate that more noise is propagated to the output than 
that would be actually propagated. However, the actual 
vectors are generally not known in situations such as 
static noise analysis (SNA) and using worst case vector en- 
sures acceptable noise in the operation of an IC for other 
vectors. Therefore, at least in SNA type situations, reduced 
computational requirements are obtained, in addition to 
using a desirable vector for noise analysis. 
[0051] In one embodiment, the parameters thus characterized 
include noise propagation (NP) data (characteristics) and 
noise immunity data/curve. The noise immunity data indi- 
cates whether an input glitch of a given height and a 
givenwould cause an output glitch exceeding a pre- 
specified threshold (e.g., 10 percent of Vdd), in which 
case the result is deemed to be a failure (generally requir- 
ing redesign of the corresponding input net) and a suc- 
cess otherwise. A noise immunity curve (NIC) may be gen- 
erated from the noise immunity data, with the NIC map- 
ping an width and height combination (of an input glitch) 
to a success or failure result, thereby enabling the analy- 
sis of an IC to be performed quickly later. Each point on 



the NIC may be referred to as a immunity transition point. 
A NIC may be generated only for the worst case vector in 
an embodiment to reduce computational requirements. 

[0052] The NP data on the other hand indicates the height and 
width of an output glitch in response to an input glitch of 
a specified height and width. A NP curve may be gener- 
ated for each width of interest. Each NP curve indicates 
the corresponding height of an output glitch in response 
to an input glitch of a given height. The NP data may also 
be generated based on only the worst case vector (for a 
given arc of a cell). As a result, the computational require- 
ments for characterization of a cell library may be re- 
duced, as noted above. 

[0053] Further reduction in computational requirements for de- 
termining NP parameters and noise immunity parameters 
may be achieved, as described below with examples. 

[0054] vvith respect to determining immunity transition point on 
a NIC, in general, a search range (with upper bound and 
lower bounds for height, while fixing width) is first de- 
fined and simulations are performed within the search 
range to predict the boundary point of the failure area and 
success area for a NIC. Various aspects of the present in- 
vention enable the search range to be narrowly defined. 



thereby reducing the number of simulations as described 
below. 

[0055] According to an aspect of the present invention, the 

search range is reduced by fitting pre-determined immu- 
nity transition points into a curve and predicting the 
height of the immunity transition point sought to be de- 
termined. The upper and lower bounds are respectively 
set to equal more and less than the predicted height by a 
specific value (e.g., 15 points of the resolution of the 
search), and the range may be searched to determine the 
exact height of the immunity transition point. Due to the 
predicted height and thus reduced search range, the 
number of simulations may be reduced. 

[0056] According to another aspect of the present invention, two 
input glitches of the same width but of heights separated 
by a step size of the NP curve/data are determined such 
that one input glitch provides a failure result and another 
glitch provides a success result. The two input glitches 
may be advantageously determined while generating NP 
data, and thus may not require additional computations. 
The search range would contain the (height corresponding 
to) immunity transition point due to the use of the NP 
data. The number of simulations may be reduced as a re- 



suit. 

[0057] According to one more aspect of the present invention, if 
multiple immunity transition points are already known, 
the points may be fitted into a curve and the height of the 
immunity transition point sought to be determined is pre- 
dicted (for a given width) by appropriate mapping on to 
the curve. The heights corresponding to the two adjacent 
(one for narrower width and another for wider width) 
points (i.e., immunity transition points already deter- 
mined) used as the bounds for the search range. As the 
search range would accurately contain the immunity tran- 
sition point sought to be determined, the number of sim- 
ulations can be reduced. 

[0058] The overlap of the search ranges noted above can be ad- 
vantageously used to further reduce the search range. 

[0059] The computational resource requirements for NP curves 
may also be reduced bythat the height of the output 
glitches for corresponding height of the input glitches 
would remain substantially same outside of a transition 
region due to either no glitch propagation or causing a 
definite transition at the output pin. search may be per- 
formed to identify the boundaries of the transition region, 
and the height corresponding to the boundary points may 



be used to set the values of the points outside of the 
transition region. Computational requirements are re- 
duced as a result. 

[0060] Another aspect of the present invention enables the start 
search point to be reliably placed in the transition region. 
Assuming a NIC is either available or plotted (using curve 
fitting), the height corresponding to the width for which 
the NIC relates to, is determined from the NIC, and the 
start point of the NP curve is set based on the determined 
height. Since NIC by definition provides (output glitch) 
height of a pre-specified percentage of the maximum 
possible height, the determined height may reliably fall in 
the transition region, thereby reducing the computational 
requirements further. 

[0061] One more aspect of the present invention reduces the 

computational requirements by interleaving the computa- 
tion of noise immunity data and NP data. In an embodi- 
ment, the NP curves corresponding to the widest and nar- 
rowest widths of input glitches are determined first. The 
immunity transition points for the two widths are then de- 
termined using the NP data to reduce the search range, as 
described above. A third NP curve may then be deter- 
mined corresponding to an intermediate (middle) width. 



The average of the heights corresponding to the two de- 
termined immunity transition points may be used as a 
start point for the third NP curve. 

[0062] The data from the third NP curve may then be used to de- 
termine the corresponding (of the same width) immunity 
transition point. The three immunity transition points may 
be used to generate an approximate NIC by using curve 
fitting techniques. The curve thus generated is then used 
to set the start points for the remaining NP curves. The NP 
curves in turn are used in generating the immunity transi- 
tion points, and thus the accurate NIC. The determination 
of NP curves and immunity transition points for the same 
width can be interleaved for increasing the accuracy of 
various predictions, and thus reducing the computational 
resource requirements. 

[0063] More aspects of the present invention reduce computa- 
tional requirements in noise immunity characterization of 
sequential cells. In general, it needs to be known whether 
input glitches of various heights and width would either 
cause a transition (a change of state) on the output pin or 
cause an output glitch to exceed a pre-specified thresh- 
old, in which case a failure result is deemed to be ob- 
tained. Otherwise, a success/pass result is deemed to be 



obtained. For each input width and height combination, 
the clocl< may need to be swept in a range (exceeding the 
width), and such sweep requires multiple simulations with 
the clock transition point being placed at various places in 
relation to the data transition time. 

[0064] An aspect of the present invention determines whether 
the sweep (and corresponding simulations) can be 
avoided by simulating each cell of interest with an input 
transition waveform having a height equaling the height 
of an input glitch of interest and providing infinite (i.e., 
sufficient time for any practical sequential element to 
latch) setup and hold times. If the result of the simulation 
indicates a success result (i.e., not sufficient noise propa- 
gated), sweep can be avoided for all widths of the input 
glitch for the same height. Thus, the computational re- 
quirements can be reduced. 

[0065] The computational requirements for characterization of 
sequential cells can be further reduced by using the vari- 
ous techniques which use other immunity transition 
points already determined, to reduce the search range of 
a height of the glitch. 

[0066] Several aspects of the invention are described below with 
reference to examples for illustration. It should be under- 



stood that numerous specific details, relationships, and 
methods are set forth to provide a full understanding of 
the invention. One skilled in the relevant art, however, will 
readily recognize that the invention can be practiced with- 
out one or more of the specific details, or with other 
methods, etc. In other instances, welLknown structures or 
operations are not shown in detail to avoid obscuring the 
invention. 
[0067] Example Circuit 

[0068] Figure lA is a block diagram of an example circuit used to 
illustrate several aspects of the present invention. Circuit 
100 is shown containing cells 110-A through 110-D and 
capacitor 130. Cell 110-A is shown connected to cell 
110-B by net 113-1, and cell 110-C is shown connected 
by net 113-2. 

[0069] Capacitor 130 represents a coupling capacitance between 
nets 113-1 and 113-2 and signal transitions on one net 
can cause a glitch on another net due to coupling through 
capacitor 130. For illustration, in the below description, it 
is assumed that 113-1 is the aggressor net and 113-2 is 
the victim net. The potential effect on the cell(s) con- 
nected to the victim net may be measured by 'noise im- 
munity' and 'noise propagation' parameters noted above. 



Various aspects of the present invention reduce the com- 
putational requirements in characterizing the parameters 
as described below in further detail. 

[0070] jjj jsjoise Immunity Parameters 

[0071 ] /// ^ General Introduction 

[0072] As noted above, noise immunity parameters generally in- 
dicate whether a glitch of a given width and/or height on 
a pin would cause the signal at an output pin of a cell to 
be altered by more than a pre-specified threshold. In gen- 
eral, a glitch at the input pin has to be greater than a cer- 
tain fixed voltage to cause a failure (i.e., response on the 
output pin exceeds a pre-specified threshold). However, a 
glitch with a large height may still not cause any failure, if 
the glitch width is very small. This is mainly because noise 
failure is related to input noise glitch energy and this en- 
ergy is proportional to the area under the glitch wave- 
form. 

[0073] In addition, the output response of a cell might be differ- 
ent for different values for the input pins of the cell. For 
example, with reference to Figure IB, cell 110-D is shown 
with input pins 113-2 and 111-1 through 111-N, and it 
may be desirable to ensure that pin 113-2 (as impacted 



by pin 113-1) is immune to noise for any combination of 
tlie N-bit values wliicli would cause the bit value on out- 
put pin 119 to change if a logical value on pin 113-2 is 
changed. Each such N-bit value, which would cause the 
bit value on output pin 119 to change, may be termed as 
a vector. Thus, the output response may be different for 
different vectors. 
[0074] Accordingly, in one prior approach, a noise immunity 
graph may be generated by applying a glitch of shape 
shown in Figure 2A, and generating the noise immunity 
curves (NIC) shown in Figure 2B for the corresponding 
vectors. With respect to Figure 2A, the glitch is shown 
with a height (voltage) H and width (time) W. Figure 2A is 
shown containing a rising glitch for illustration and the 
description is provided substantially with respect to rising 
glitches only for conciseness. However, similar characteri- 
zation is generally performed with reference to falling 
glitches as well. The height and width of a glitch can be 
varied to generate the curves of Figure 2B as described 
below. 

[0075] With reference to Figure 2B, the X-axis represents width 
(W) and Y-axis the height (H). In an embodiment, the Y- 
axis contains a range from 0 to twice the Vdd value. 



wherein Vdd represents the supply voltage. The range of 
values of interest on the X-axis may be determined based 
on typical values expected from design/implementation 
considerations, and are provided as inputs for characteri- 
zation. Curves 251-254 represent the immunity curves 
corresponding to four vectors. 

[0076] Each immunity curve separates the points of failure (i.e., 
the input glitches which would cause the output glitch to 
exceed a pre-specified threshold level) and the points 
where there are no failures. The points containing all the 
failure points may be viewed as being contained in a 'fail- 
ure area", and the points where there is no failure may be 
viewed as being contained in a 'safe/success area'. The 
points in the boundary between the two areas may be re- 
ferred to as immunity transition points. Thus, each immu- 
nity curve is defined by the corresponding immunity tran- 
sition points. With reference to curve 251 of Figure 2C, all 
the points above curve 251 represent failure area, all the 
points below curve 251 represent success area and all the 
points on curve 251 represent immunity transition points. 

[0077] For illustration, it is assumed that the cell sought to be 
characterized represents full adder 300, having three in- 
puts (301-303) and two outputs (sum 391 and carry 392) 



as shown in Figure 3. It is assumed that an arc between 
input pin 301 and output pin 391 (carrying the sum re- 
sult) is considered for characterization. For an arc be- 
tween inputs and sum 391, four vectors would present, 
which correspond to the four values 00, 01, 10, and 11 on 
pins 302 and 303 respectively. Curves 251-254 represent 
the immunity curves corresponding to the four vectors. 

[0078] Thus, in a prior approach, the immunity transition points 
on each curve are determined by first fixing an input vec- 
tor and one of height and width, and varying the other 
one of height/width using binary search techniques until 
the boundary of safe area and failure area is determined 
with a desired accuracy. If the height of the output glitch 
exceeds the pre-specified threshold, a 'failure result' is 
said to be obtained, and a 'success result' is said to be 
obtained otherwise. Thus, the immunity transition point is 
found when adjacent input glitches differing by the de- 
sired accuracy produce different results. 

[0079] Each point on the NIC may be thus determined for a given 
vector. The technique may be repeated for each curve. 
Once the curves are determined, analysis can be quickly 
performed to determine any potential failures in chip level 
noise analysis. It may be appreciated that substantial 



computational resources may be consumed in determin- 
ing the points on all the curves. 

[0080] For example, assuming that the maximum input glitch 
height is 2*Vdd and NIC needs to be measured to an ac- 
curacy of 0. 1% of Vdd, then total number of possible 
points in the binary search range (of 0 to 2*Vdd) would be 
2000 and hence the number of binary search iterations 
(simulations) required equals 10. Assuming NIC is per- 
formed for five vectors, eight loads (fan-out), eight 
widths, three arcs, two glitch types (falling and rising), 
three PTV (process, temperature, voltage) conditions and 
10 binary search iterations, the number of computations 
(Nnic) required to determine NIC for a cell equals (5 * 8* 8 
* 3* 2* 3*10 = 57600) in one prior approach, wherein 
'*'represents multiplication operation. 

[0081] A library may contain 1000 cells, and thus the total num- 
ber of simulations is of the order of 50 millions, which re- 
quires substantial computational resources. Various as- 
pects of the present invention reduce such computational 
requirements as described below. 

[0082] ni.B. Method 

[0083] Figure 4 is a flow chart illustrating the manner in which 
computations may be reduced while characterizing noise 



immunity parameter according to an aspect of tlie present 
invention. Tlie flow chart is described with reference to 
Figure 5 (having width and height of glitches on X and Y 
axis respectively) for illustration. The flow-chart begins in 
step 401, in which control immediately passes to step 
410. 

[0084] In step 410, the input vectors of interest for an input and 
output path combination of a cell being characterized are 
determined. In an embodiment, a combination of bit val- 
ues for the other pins is considered an input vector only if 
a change in logic value on the input pin would cause a 
change in logic value on the output pin of the cell. The in- 
put vectors can be determined by analyzing the truth table 
of cell to be characterized. For illustration, it is assumed 
that three vectors are identified. 

[0085] In step 430, a small number of widths or heights (one di- 
mension/coordinate) of the glitch are selected. With re- 
spect to the graph of Figure 2B, such a step may be im- 
plemented by selecting a small number of points on either 
X-axis or Y-axis. In an embodiment, two points are se- 
lected towards the respective ends and one point is se- 
lected towards the middle, where the curves are likely to 
inflect/curve substantially. For example, with reference to 



Figure 5, points 510, 512 and 511 respectively represents 
tlie widtlis at lower end, higher end and middle where the 
inflection points are expected for the immunity curves. 
[0086] In step 450, the height of immunity transition points cor- 
responding to the three widths are determined for each 
immunity vector. The height (or other dimension for the 
immunity transition point, in general) can be determined 
using various well-known techniques. In one embodiment, 
a binary search technique is employed to determine the 
height for each point. The number of searches/simula- 
tions can be reduced further according to various aspects 
of the present invention, as described in further detail in 
sections below. 

[0087] With respect to Figure 5, points 531, 532 and 533 are de- 
termined respectively for the three vectors for first width 
510, points 541, 542 and 543 are determined respectively 
for the three vectors for third width 511, and points 551, 
552 and 553 are determined respectively for the three 
vectors for second width 512. 

[0088] In step 460, curve (representing approximate noise im- 
munity curve) corresponding to each vector is plotted 
based on the immunity transition points. Curve fitting 
techniques may be employed to generate each noise im- 



munity curve. In one embodiment, the curves are gener- 
ated using Levenberg-IVIarquarndt teclinique described, 
for example, in a text book entitled, "Numerical Recipes: 
The Art of Scientific Computing" by B.P. Flannery, S. A. 
Teukolsky and W. T. Vetterling, (ISBN_0_521_43108_5), 
Cambridge University Press. An example of such an ap- 
proach is described in sections below. 

[0089] With reference to Figure 5, immunity curves 521, 522 and 
523 are respectively shown plotted using respective sets 
of points {521, 531, and 541}, {522, 532, and 542} and 
{523, 533, and 543}. It should be understood that the by 
including more points in steps 430 ad 450, more accurate 
curves can be generated. However, the computational re- 
quirements are correspondingly enhanced. 

[0090] In step 470, the vector corresponding to the lowest noise 
immunity curve (i.e., with least height for each given 
width) is set as the worst case vector. With reference to 
Figure 5, vector 521 at the lowest position among the 
three vectors is set as the worst case vector. 

[0091] In step 480, immunity transition points for the remaining 
widths corresponding to only the worst case vector are 
determined. That is for each point in one dimension, the 
value in the other dimension is determined. The corre- 



spending noise immunity points are plotted to generate 
the noise immunity curve accurately. 

[0092] In step 490, an IC is analyzed for cross-talk noise using 
the immunity curve corresponding to the worst case vec- 
tor determined in step 480. In general, the glitch (on an 
input pin) width and height would be determined/com- 
puted during the noise analysis, and the worst case im- 
munity curve is examined to determine whether the com- 
puted output glitch is in the failure area or safe area. If 
there is a failure result (or result falls in the failure area), 
the portion of the circuit affecting the signal on the corre- 
sponding net would need to be redesigned. 

[0093] In an embodiment, the same noise immunity curve is used 
to determine failures/violations, irrespective of the spe- 
cific vector applied to the inputs. If a cell is deemed to be 
acceptable in the analysis with respect to the worst case 
vector, the cell may be deemed to be acceptable with 
other vectors as well. However, as noted above, the actual 
vectors are generally not known in situations such as 
static noise analysis (SNA) and using worst case vector en- 
sures acceptable noise in the operation of an IC for other 
vectors. The method ends in step 499. 

[0094] The approach of above was described assuming that the 



noise immunity curves do not cross one anotlier, as 
shown in Figure 5. However, tlie noise immunity curves 
can cross over and tlie approach described above can be 
extended as described below. 

[0095] Noise Immunity Curves with Cross-overs 

[0096] Figure 6 depicts three noise immunity curves 601, 602 
and 603, with the curves 601 and 602 crossing. Curve 
603 represents the worst-case immunity curve (and is not 
crossing other curves). In such a case (in which the worst 
case curve is not crossing), the approach of Figures 4 and 
5 can be used. The description is continued with respect 
to a scenario there is a crossing of curves at the lowest 
level, and thus there is not one curve which can be termed 
as a worst case curve. 

[0097] Figure 7 depicts a scenario in which two curves (702 and 
703) at the lowest level are crossing over at time point 
713. In such a case, in one embodiment, each of the two 
curves is treated as a worst case curve in all further analy- 
sis. That is, characterization data may be generated for 
both curves. 

[0098] In an alternative embodiment, the width until which the 
two curves do not cross over, the lower curve 703 is used 
as the worst case curve and after the cross-over, two 



curves 702 and 703 are used as the worst case curves. 
Due to the use of such worst case transition points from 
both the curves, the approach may reliably indicate at 
least problems which are likely to occur. 

[0099] Noise Propagation (NP) Parameters 

[0100] /v: A. General Introduction 

[0101] As noted above, noise propagation parameters specify 

whether an output pin of a cell would switch (and possibly 
cause a wrong value to be latched in a subsequent cell) in 
response to a glitch of a given width and height. Noise 
propagation is measured for combinational cells. In gen- 
eral, NP data is measured by injecting an input glitch of 
l<nown width and l<nown height and measuring corre- 
sponding output glitch height and output glitch width. In 
one prior approach, the NP data is plotted as described 
below with reference to Figures 8A and 8B. 

[0102] Figure 8A is a graph illustrating the input glitches applied 
to a buffer and Figure SB depicts the corresponding re- 
sponses at an output pin. All input glitches 810-1 
through 810-8 are shown having the same width and the 
output responses of Figure 8B correspond to a single vec- 
tor. 



[0103] As niay be appreciated, up to some threshold value of in- 
put voltage, the buffer blocks the input from being propa- 
gated to the output. In other words, the buffer does not 
respond to certain threshold heights in case width is 
fixed, as shown in Figure 8B. Accordingly,corresponding 
to input glitches 810-1 through 810-3 are not shown in 
Figure 8B. On the other hand, output responses 820-4 
through 820-8 are generated in response to input glitches 
810-4 through 810-8. 

[0104] Further, the input peaks for glitches 810-1 through 

810-8 are respectively shown at 830-1 through 830-8, 
and the output peaks are shown as points 840-4 through 
840-8 when the buffer passes the input glitch in the case 
of output responses 820-4 through 820-8. 

[0105] Figure 9 contains a graph with the peak values of input 

glitches shown on X-axis and the corresponding peaks of 
the output responses shown on Y-axis. The peak value is 
shown as zero until the input peak equals or exceeds 
810-4 (as shown at point (830-4, 840-4)), and saturates 
at voltage level Vdd corresponding to input peak 830-7 
(as shown at point (830-7, 840-7)). Even if the input peak 
is increased after that point, the output peak remains the 
same, as shown by the flat portion. Similarly, the flat por- 



tion below represents a situation in which the input glitch 
is blocked until the peak exceeds certain value (911). 

[0106] Voltage level 917 represents the threshold voltage for the 
noise immunity analysis, and represents 10% of the supply 
voltage (Vdd). Input peak value on X-axis corresponding 
to voltage level 917 represents immunity transition point, 
which may be used to determine NIC as described in sec- 
tions below with reference to Figure 20. 

[0107] It should be appreciated that the curve of Figure 9 may be 
referred to as a noise propagation (NP) curve, and is only 
for a combination of one width and one vector. In one 
prior embodiment, a similar NP curve is generated for var- 
ious combinations of widths and vectors before analysis 
for cross talk noise is performed for an integrated circuit 
(IC). 

[0108] During analysis, the width and height of an input glitch at 
an input pin/net may be provided by the analysis tool, and 
the NP curve (data) may be analyzed to determine the ex- 
tent to which the input glitch is propagated to an output 
pin. The glitch that would be presented on output pin 
could be used in analyzing the behavior of subsequent 
cells. From the above, it may be appreciated that substan- 
tial computational resources may be required to deter- 



mine the curves (as in Figure 9) corresponding to various 
combinations of glitcli width and vector values. 

[0109] For example, assuming that the maximum input glitch 

height is 2*Vdd and NP step size equals 10% of Vdd, each 
NP curve would contain 20 points. Assuming NP chracteri- 
zation is performed for five vectors, eight loads (fan-out), 
eight widths, three arcs, two glitch types (falling and ris- 
ing), three PTV (process, temperature, voltage) conditions 
and 20 points in each NP curve, the number of computa- 
tions (Nnp) required to chracterize NP for a cell equals (5 * 
8* 8 * 3* 2* 3*20 = 115200) in one prior approach, 
wherein '*'represents multiplication operation. 

[0110] A library may contain 1000 cells, and thus the total num- 
ber of simulations is of the order of 100 millions, which 
requires substantial computational resources. An aspect 
of the present invention reduces the computational re- 
quirements as described below with reference to Figure 
10. 

[0111] lY B. Reducing Computational Requirements for Noise Propagation 
Analysis 

[0112] Figure 10 is a flow chart illustrating the manner in which 
computations may be reduced while performing noise 
propagation analysis according to an aspect of the present 



invention. The flow-chart begins in step 1001, in which 
control immediately passes to step 1010. 
in step 1010, the worst case vector for an input pin and 
output pin combination, which would propagate most 
noise from the input pin to the output pin among all the 
vectors of interest, is determined. Such a worst case vec- 
tor may be determined using various techniques. 
[0^^^] In one embodiment, the worst case vector set in step 470 
is determined as the worst case vector for step 1010 as 
well. Such an approach is based on the fact that there is 
commonality between noise propagation analysis and 
noise immunity analysis given that a cell generally gets 
configured (e.g., transistors being off/on) based on the 
input vector, and the signal propagation occurs from sim- 
ilar configuration in the case of both the analysis. In addi- 
tion, since by definition, failure result for a NIC is placed 
at specific points (corresponding to 10% of maximum 
height of the output glitch) on the NP curves, such use of 
the same vector as a common worst case vector may be 
justified. 

[0115] In step 1050, the NP data curve (similar to in Figure 9) is 
determined for the worst case vector of step 1010 for 
each (glitch) width of interest. Each point of each NP data 



curve may be determined using a known way. 

[01 16] In step 1080, an integrated circuit (IC) is analyzed using 
tlie NP data curve corresponding to the worst case vector. 
That is, once a input glitch of a given width and height is 
known, the NP data curve corresponding to the width may 
be selected first, and the height (X-axis) of the input 
glitch may be used to determine the expected height of 
the output glitch from the NP noise data. The expected 
height of the output glitch can be used to analyze the IC 
further. The flow chart ends in step 1099. 

[01 1^] It may be appreciated that the computational complexity 
is substantially reduced since NP data (curve) is deter- 
mined for only one vector. In one scenario, the number of 
simulations were reduced by more than half compared to 
a scenario in which parameter data related to all the vec- 
tors is determined. In addition, as the worst case vector is 
determined in common for both immunity and propaga- 
tion analysis, the additional overhead is reduced. As noted 
above, at least in static analysis type situations, the data 
corresponding to worst case vectors may be acceptable 
(or even desirable). 

[0118] y Reducing Computations in Determining Immunity Transition 
Points 



[0119] y ^ General Introduction 

[0120] As described above with reference to Figure 2A, immunity 
transition points form an immunity curve separating tlie 
point of failure ("failure area") and the "safe area". There is 
often a need to determine the immunity transition points. 
For example, as noted above, with reference to step 480, 
immunity transition points may need to be determined for 
the worst case vector with the lowest noise immunity. 

[0121] In one prior embodiment, a 'search' is performed between 
0 and Vdd voltage (maximum possible height for input 
glitches) for each immunity transition point. In general, a 
search entails multiple simulations until the immunity 
transition point is determined with a desired accuracy 
(i.e., a change of the immunity transition point by the de- 
sired accuracy should change the result of the noise im- 
munity analysis). Simulations require computational re- 
sources, and it is desirable to reduce the number of simu- 
lations. An aspect of the present invention enables such 
reduction in the search range (and thus the number of 
simulations), as described below in further detail. 

[0 1 22] Determining Noise Immunity Curve 

[0123] Figure 11 is a flow chart illustrating the manner in which a 



noise immunity curve may be determined according to an 
aspect of tlie present invention. Tlie flow cliart is de- 
scribed with reference to the graph of Figure 12, which 
contains the width of an input glitch on X-axis and the 
height on the Y-axis. The flow chart begins in step 1101, 
in which control immediately passes to step 1110. 

[0124] In step 1110, a variable T representing the width number 
presently being determined, is set to 1 (indicating that the 
first width is being computed). Variable A, representing 
the desired accuracy is set to (Vdd * 0.0001), wherein Vdd 
represents the supply voltage, the lower bound (LB) is set 
to 0, and the upper bound (UB) is set to (2 * Vdd). 

[0125] In step 1120, a determination is made as to whether i is 
greater than 3. Control passes to step 1130 if i greater 
than 3, and to step 1150 otherwise. As will be appreci- 
ated, after the first three iterations (corresponding to i 
values of 1, 2 and 3), the initial search range for subse- 
quent widths is determined according to an aspect of the 
present invention to minimise the computational require- 
ments. 

[0126] The steps of 1150, 1160 and 1170 operate to perform a 
search in the initial range defined variables UB and LB for 
each value of i.step 1150, the output response corre- 



spending to the input glitcli of widtli W[i] is examined. 
Tlie output response is generated by simulation using 
tools such as SPICE simulators available in the market 
place. 

[0127] In step 1160, a determination is made as to whether the 
immunity transition point is found with accuracy A. An 
immunity transition point would be deemed to be found if 
a change of value by amount A would change the result of 
analysis from failure to pass or vice versa. Control passes 
to step 1180 if immunity transition point is found, and to 
step 1170 otherwise. 

[0128] In step 1170, the UB and LB values are adjusted according 
to the results of the simulation. Assuming a binary search 
technique, depending the failure or pass (not failure) re- 
sult of the prior simulations, the UB and/or LB values are 
adjusted to correspond to the binary search principle, and 
may be performed in a known way. Control again passes 
to step 1150 with the resulting UB and LB values. 

[0129] Steps 1130 and 1140 operate to reduce the computational 
requirements after the first three data points of a immu- 
nity curve are determined. As may be appreciated from 
the description above, the search for each of the first 
three noise immunity points is in the range 0 and (2*Vdd) 



due to the initialization in the step 1110. However, due to 
the operation of steps 1130 and 1140 the search range is 
reduced to a narrower range, thereby reducing the com- 
putation requirements, as described below. 

[0130] For illustration, with reference to Figure 12, it is assumed 
that points 1201, 1202 and 1203 are determined in the 
first three iterations. The manner in which the remaining 
immunity transition points are determined according to 
various aspects of the present invention, is described be- 
low in further detail. 

[0131] Continuing with reference to Figure 11, in step 1130, a 
variable Ynic is measured by fitting into an immunity 
curve the heights measured thus far (only three when 
control passes to step 1130 the very first time), and by 
mapping the width of the present iteration into curve. 
Such an approach provides a fairly good prediction of the 
height for noise immunity analysis purpose since it has 
been observed that noise immunity curves approximately 
follow a hyperbolic equation for each vector. Various 
techniques can be used for curve fitting and Levenberg- 
Marquardt technique noted above is used for the curve 
fitting. 

[0132] vvith reference to Figure 12, it is assumed that curve 1220 



is generated by appropriate curve fitting technique and to 
determine tlie glitcli lieight corresponding widtli W2, vari- 
able Ynic is set to the Y-coordinate for the point 1204. 
The search range is set between Y-coordinates 1231 and 
1232 for the glitch width of W2 as described below. 

[0133] In step 1140, the search range is set to a narrow band 

around the computed Ynic. For illustration, it is assumed 
that the search range is set between (Ynic+ 15*A) and 
(Ynic - 15*A) by setting the variables UB and LB appropri- 
ately. Control then passes to step 1150. As noted above, 
steps 1150, 1160 and 1170 operate to determine the 
minimum height at which the output signal of a cell ex- 
ceeds a pre-specified threshold (for the present width). 

[0^34] In step 1180, variable i is incremented by one to process 
the next width. In step 1190, the value of i is compared 
with N (total number of widths) and control passes to step 
1199 if i is greater than N, and to step 1110 otherwise. 
The method ends in step 1199. It should be appreciated 
that the curve can be refined further as each immunity 
transition point is accurately determined, thereby enhanc- 
ing the accuracy of the search range in determining the 
next immunity transition points. 

[0135] Due to the use of a smaller search range in steps 1130 



and 1140, the computational requirements in generating 
a noise immunity curve may be reduced. For example, 
without using steps 1130 and 1140, the total number of 
possible points to search will be 2000 (binary search 
range 0 to 2VDD), and 10 binary search iterations will be 
required to get NIC data for a given input glitch. On the 
other hand, by using steps 1130 and 1140, the number of 
iterations can be reduced to 5 (since 30 points need to be 
searched), thereby reducing the computational require- 
ments by half. Various aspects of the present invention 
further reduce the computational requirements, as de- 
scribed below. 

[0136] yj^ Improving Accuracy of Search Range to Determine Immunity 

Transition Points 
[0137] VT. A. General Introduction 

[0138] One problem with the approach of Figures 11 and 12 de- 
scribed above is that the prediction of possible value/ 
height of a immunity transition point for a given width 
(and thus the range), is based on mathematical modeling 
of the expected NIC. The range determined accordingly 
may not contain the actual immunity transition point, and 
error recovery approaches may need to be implemented to 
search outside of the search range determined by steps 



1130 and 1140. The computational requirements are en- 
hanced in such scenario, and it is generally desirable to 
determine a search range accurately. 

[0139] An aspect of the present invention takes advantage of 

various pieces of other information to improve the accu- 
racy of search range as described below in further detail. 
First, the manner in which NP data can be used to improve 
the search range accuracy, is described below. 

[0140] Principle: Using NP Data 

[0141] Figure 13A is a timing diagram of input glitches and cor- 
responding output glitches while determining NP data. 
The timing diagram illustrates the manner in which the 
search range can be determined accurately for a given 
width of input glitches. The diagram is shown containing 
input glitches 1301-1304 of respective heights 
1311-1314, and the corresponding output responses 
1331 - 1334. For illustration, it is assumed that the 
height of output response 1333 exceeds the pre-specified 
threshold level (10% of Vdd) and the height of output re- 
sponse 1332 is below the pre-specified threshold level. 

[0142] According to an aspect of the present invention, the 

search range is set between 1312 and 1333. In an em- 
bodiment, the heights of input glitches are spaced apart 



by 10% of Vdd, and thus the search range equals 10% of 
Vdd. The search range can be determined within a few 
simulations, for example, by using the approach de- 
scribed below with reference. However, the range would 
be accurate given that by definition, the height sought to 
be determined by the noise immunity analysis falls within 
this range. Therefore, the search range can be determined 
accurately by using this principle. A method using the 
principle is described below. 
[01 43] vr.C Method: Using NP Data 

[0144] Figure 13B is a flowchart illustrating the manner in which 
the search range to determine an immunity transition 
point can be determined accurately according to an aspect 
of the present invention. The method begins in step 1351, 
in which control passes to step 1355. 

[0145] In step 1355, two input glitches of the same/given width 
having a difference of height equaling the step size of the 
NP curve of the same width, with one input glitch causing 
failure result and the other input glitch causing success 
result in noise immunity analysis, are determined. The two 
input glitches may correspond to glitches 1302 and 1303 
described above since the two input glitches have heights 
differing by the step size of the corresponding NP curve. 



[0146] In step 1365, the two bounds of the search range are set 
to equal the respective height of the two input glitches. 
Thus, in the example of above, the upper and lower 
bounds of the search range equals value corresponding to 
points 1313 and 1312 respectively. Techniques such as 
binary searches may be applied to determine the precise/ 
accurate immunity transition point as described above 
with reference to steps 1150, 1160 and 1170 of Figure 
11. As also noted there, the range determination of Figure 
13B and search of Figure 11 can be repeated for different 
widths of interest. The method ends in step 1399. 

[0147] Thus, the approach of Figure 13B enables accurate deter- 
mination of a search range. Other approaches can also be 
used to further constrain the search range as described 
below. 

[01 48] yj£)_ Using Data from Determined Immunity Transition Points 

[0149] Figure 14A is a flowchart illustrating the manner in which 
data related to determined immunity transition points may 
be used to determine the search range accurately while 
generating a noise immunity curve. The approach is based 
on a recognition that larger (or at least not less than) 
height is required to cause a failure result when a nar- 
rower glitch is used. The flowchart is described with refer- 



ence to Figure 14B (containing time/width on X-axis and 
lieiglit of input glitcli on Y-axis) for illustration. The 
method of Figure 14A begins in step 1401, in which con- 
trol immediately passes to step 1405. 
[0150] In step 1405, a first immunity transition point is deter- 
mined, with the first immunity transition point having a 
width greater than the width of the immunity transition 
points sought to be determined. Similarly, in step 1415, a 
second immunity transition point is determined, with the 
second immunity transition point having a width less than 
the width of the immunity transition points sought to be 
determined. 

[0151] The first and second immunity transition points may be 
determined based on pre-determined information. In an 
embodiment, immunity transition points corresponding to 
narrow width and widest width are determined first by 
simulation and the immunity transition point correspond- 
ing to middle width is sought to be determined. In such a 
situation, the first and second immunity transition points 
would equal the immunity transition points corresponding 
to the narrow and widest widths. The first and second 
transition points corresponding to remaining widths may 
be determined based on the pre-determined immunity 



transition points of closest widths. 

[0152] For example, with reference to Figure 14B, assuming that 
pointsl401 and 1402 are already determined correspond- 
ing to widths 1421 and 1422 respectively, and width 1411 
represents the width for which the immunity transition 
point is sought to be determined, points 1401 and 1402 
respectively correspond to the second and first immunity 
transition point respectively. 

[0153] In step 1420, the lower bound of the search range is set 
to equal the height of the first immunity transition point, 
and in step 1425 the upper bound of the search range is 
set to equal the height of the second immunity transition 
point. Thus, the search range is between coordinate val- 
ues represented by 1431 and 1432. The method ends in 
step 1449. 

[0154] As may be appreciated, the search range according to Fig- 
ures 14A and 14B becomes more accurate as more immu- 
nity transition points are determined. The approaches of 
above (Figures 11, 13B and 14A) can be combined to de- 
termine a narrow and accurate search range as described 
below with reference to Figures 15 and 16. 

[0 1 55] yj£ Accurate and Narrow Search Range for Noise Immunity Anal- 
ysis 



[0156] Figure 15 is a flow chart illustrating the manner in which 
various approaches described above can be combined to 
generate a narrow and accurate search range. The 
flowchart is described with reference to Figures 16A and 
16B. The flowchart begins in step 1501 in which control 
immediately passes to step 1510. 

[0157] In step 1510, a first search range is determined using 
noise propagation (NP) data. The approaches described 
above with respect to Figures 13A and 13B may be used 
for such a determination. 

[0158] In step 1520, a second search range is determined based 
on the immunity transition points corresponding to adja- 
cent widths, which are higher and lower than the width for 
which the immunity transition point is sought to be deter- 
mined. The approaches described above with respect to 
Figures 14A and 14B may be used for such a determina- 
tion. 

[0159] In step 1530, a third search range is determined as the 
intersection of the first search range and the second 
search range. It may be appreciated that the third search 
range represents an accurate search range since both the 
first search range and the second search range would 
contain the immunity transition point sought to be deter- 



mined. Further refinement of tlie searcli range may be 
possible, as described below. 

[0160] In step 1540, a fourth search range is determined from 
the curve fitted using immunity transition points deter- 
mined already. The determination may be performed as 
described above with respect to Figures 11 and 12. 

[0161] In step 1560, a determination is made as to whether there 
is a overlap of the third search range and the fourth 
search range. Control passes to step 1570 if there is a 
overlap, and to step 1580 otherwise. 

[0162] Step 1570 corresponds to a situation depicted in Figure 
16A, in which the entire potential search range from 0 to 
2 * Vdd is shown as 1601. The first, second, third and 
fourth search ranges are shown as 1602, 1603, 1604 and 
1605 respectively. As can be readily appreciated, there is 
a overlap between ranges 1604 and 1605, and the overlap 
is shown as range 1606. In step 1570, range 1606 is set 
as the desired range in which further search is to be con- 
ducted to determine the immunity transition point. 

[0163] Step 1580 corresponds to the situation in Figure 16B, in 
which the entire potential search range from 0 to 2 * Vdd 
is shown as 1651. The first, second, third and fourth 
search ranges are shown as 1652, 1653, 1654 and 1655 



respectively. As can be readily appreciated, there is no 
overlap between ranges 1654 and 1655. In step 1570, 
third range 1654 is set as the desired range in which fur- 
ther search is to be conducted to determine the immunity 
transition point. The method ends in step 1599. 

[0164] It niay be appreciated that the search may be determined 
accurately and narrowly by using some of the techniques 
described above while performing noise immunity charac- 
terization. The manner in which computational require- 
ments can be reduced for NP characterization also, is de- 
scribed below in further detail. 

[0 1 65] yji^ Reducing Computations Further in Determining NP Data 

[0 1 66] yu A^ General Introduction 

[O^^y] It is helpful to note that NP data indicates the height and 
width of an output glitch in response to an input glitch of 
a given height and width. Of interest is NP curve 
1710depicted in Figure 17. NP curve 1710 represents the 
height of the output glitch (on Y-axis) for a corresponding 
height of the input glitch (on X-axis). NP curves may be 
generated for each width of interest. 

[0168] In one prior embodiment, a simulation is performed for 

each height of interest of input glitch to generate the cor- 



responding points on NP curve 1710. For example, as- 
suming tliere are 20 input peal<s of interest (i.e., number 
of points on eacli NP curve) corresponding to a step size 
of 10% of Vdd (maximum input glitch height of 2 * Vdd), a 
corresponding number of simulations are performed for 
each width (NP curve). An aspect of the present invention 
reduces the number of such simulations, as described be- 
low with respect to Figure 18. 
[0 1 69] yjiB. Reducing Number of Simulations 

[0170] Figure 18 is a flowchart illustrating the manner in which 
the number of simulations can be reduced while generat- 
ing a NP curve (for a given width) according to an aspect 
of the present invention. The flowchart is described with 
reference to Figure 17. The flowchart starts in step 1801 
in which control immediately passes to step 1810. 

[0171] In step 1810, a start point of height of input glitch (of a 
given width) is predicted. As described below with refer- 
ence to Figure 19 below, the start point should preferably 
fall between points 1701 and 1702, in which the curve is 
operating in a transition region (i.e., the output peak 
changes substantially with a change in the input peak). 

[0172] In step 1815, the predicted start point is set as a present 
height. In step 1820, the NP data (the height and output 



of the output glitch) corresponding to the present height 
is measured. In an embodiment, a simulation tool pro- 
vides the area and height of the output signal. The output 
glitch is approximated to a triangle signal, and the width 
of the output glitch is computed based on the area and 
the height. 

[01^3] In step 1830 a determination is made as to whether the 
height of output glitch is less than Vdd. If the height of 
output glitch is less than Vdd, control passes to step 
1840, otherwise to step 1850. In step 1840, the present 
height is incremented by a small value. In an embodiment, 
the present height is incremented by 0.1 volts (10% of 
Vdd). Control then transfers to step 1820. 

[0M4] Control reaches step 1850 when the input height is 

greater than or equal to the coordinate value correspond- 
ing to point 1702. In step 1850, the height of output 
glitch in NP data is set to Vdd for all remaining higher 
heights (with X-coordinate value greater than 1702) of the 
input glitch. Such an approach may be appreciated by 
noting that the output glitch height remains substantially 
same in response to increments after point 1702. In addi- 
tion, the width of the output glitches may not be needed 
in such a region. 



In step 1855, the present height is decremented by a 
small value (say 0.1 Vdd). In step 1860, the NP data corre- 
sponding to the present height is measured. In step 1875, 
the height of the output glitch is compared with zero (or 
substantially close to 0). If the height of the output glitch 
is greater than the compared value, control passes to step 
1855, otherwise to step 1885. 
[0176] Control reaches step 1885 when the input height (height 
of the input glitch) is less than or equal to the coordinate 
value corresponding to point 1701. In step 1885, the 
height of output glitch in NP data is set to zero for re- 
maining lower heights (i.e. with X-coordinate value less 
than 1701) of the input glitch. The flow chart ends in step 
1899. 

[0177] As may be appreciated, due to the operation of steps 

1850 and 1855, the number of simulations, and thus the 
computational requirements, are reduced. The description 
is continued with respect to the manner in which the start 
point noted in step 1810 can be determined. As noted 
above with respect to step 1810, the start point needs to 
be preferably placed in the region between points 1701 
and 1702. 

[0178] However, the distance between 1701 and 1702 is narrow 



and challenges may be presented to computationally pre- 
dict the start point within the region. In an embodiment, 
the transition region equals 100-200mV for a Vdd of 
1.5V. Various aspects of the present invention enable such 
a goal to be achieved with reduced computations. The 
technique may differ for the first three NP curves and the 
later NP curves, as described below. 
[0 1 79] yu c. Determining Start Point for First Three NP Curves 

[0180] Figure 19 is a flowchart illustrating the manner in which 
the start point may be placed in the transition region for 
the first three NP curves according to an aspect of the 
present invention. The flowchart is described with refer- 
ence to Figure 20. The flowchart starts in step 1901 in 
which control immediately passes to step 1910. 

[0181] In step 1910, the start point (of height of input glitch) is 
set closer to Vdd for a NP curve (2001 of Figure 20) corre- 
sponding to a narrowest width of the input glitch. In an 
embodiment, the start point corresponding to the narrow- 
est width is set equal to 0.7Vdd. 

[0182] In step 1920, the start point of height of input glitch is set 
farther Vdd (closer to Vss) for the NP curve (2003 of Fig- 
ure 20) corresponding to the widest width of the input 
glitch. In an embodiment, the start point corresponding to 



the widest width is set equal to 0.3Vdd. 
[0183] The reason underlying such selection of the start points in 
steps 1910 and 1920 may be appreciated by noting that a 
failure result requires certain energy level (proportionate 
to area, which can be computed by multiplying of height 
and width) of the input glitch, and thus a higher height is 
chosen for narrow width and lower height is chosen for 
wider widths. 

[0184] In steps 1930 and 1940, the NP curves corresponding to 
the narrower width and widest width are respectively de- 
termined. The determination may be performed as de- 
scribed above with reference to Figure 18. NP curves 2001 
and 2003 are respectively generated as a result. 

[0185] In step 1950, a first value is set equal to the input peak 
causing a failure result (in noise immunity characteriza- 
tion) from NP curve corresponding to the narrower width. 
With reference to Figure 20, the first value equals the in- 
put peak at point 2011, corresponding to an output peak 
equaling the threshold value of (0.1 * Vdd). 

[0186] In step 1960, a second value is set equal to input peak 
causing a failure result from NP curve corresponding to 
the wider width. The second value equals the input peak 
at point 2013 in Figure 20. It may be noted that the input 



peaks of steps 1950 and 1960 may be determined accu- 
rately for noise immunity analysis, and such accurate val- 
ues may be advantageously used, as will be clearer from 
the description below with reference to Figure 22. 

[0187] In step 1970, the start point for the third NP curve 2002 

(corresponding to the middle width) is set equal to the av- 
erage of the first value and the second value. In step 
1980, the NP curve is determined for the middle width, as 
described above with reference to Figure 18. The method 
then ends in step 1999. The description is continued with 
reference to the manner in which the NP curves for other 
widths may be determined according to an aspect of the 
present invention. 

[0 1 88] yn.D. Determining Start Point for the Remaining NP Curves 

[0189] Figure 21 is a flowchart illustrating the manner in which 
the start point may be placed in the transition region for 
the remaining (after first three) NP curves according to an 
aspect of the present invention. The flowchart is described 
with reference to Figure 12. The flowchart starts in step 
2101 in which control immediately passes to step 2110. 

[0190] In step 2110, a noise immunity curve is plotted for a given 
vector. Some of the immunity transition point may be de- 
termined accurately as described in sections above and 



the other portion of the curve may be generated by the 
curve fitting techniques noted above, 
in step 2120, the start point is predicted based on the 
height on the NIC corresponding to the width for which 
the NP curve is being generated. For example, assuming 
that a NP curve is being determined between NP curves 
2001 and 2002, that the corresponding immunity transi- 
tion points respectively equal 1201 and 1202, and that 
the NP curve sought to be determined corresponds to a 
width (X-axis of Figure 12) of immunity transition point 
1204, the start point would be set equal to Y-coordinate 
of point 1204. 

[0192] In step 2150, the NP curve is determined corresponding to 
the given width using the approaches described above 
with respect to Figure 18. The method ends in step 2199. 

[0193] The description is continued with reference to a broad ap- 
proach which takes advantage of various techniques de- 
scribed above to further reduce the computational re- 
quirements. 

[0 1 94] vht; Determining NP Curves and NIC Together 

[0195] Several aspects of the present invention reduce computa- 
tional requirements in characterizing NP data by using 
data available from the NICs, and also in characterizing 



NICs by using NP data as described now. Figure 22 is a 
flowcliart illustrating the manner in which the NP curves 
and the NICs may be determined together according to an 
aspect of the present invention. The flowchart starts in 
step 2201 in which control immediately passes to step 
2210. 

[0196] In step 2210, the NP curves corresponding to the narrow- 
est width and widest width for all vectors are determined. 
The approaches described above with reference to Figures 
18, 19 and 20 may be used to determine the curves for 
each vector. 

[019^] In step 2215, the accurate immunity transition points cor- 
responding to the narrowest width and widest width may 
be determined for all vectors based on NP curves deter- 
mined in step 2210. The approach described above with 
respect to Figure 13B may be used to determine a narrow 
search range using the data available from step 2210, and 
the approach described with reference to Figure 11 may 
be used to determine the specific accurate immunity tran- 
sition point within the search range. Due to the use of the 
information available from step 2210, the accurate immu- 
nity transition points may be determined and computa- 
tional requirements are reduced in step 2215. 



[0198] In step 2220, the NP curve corresponding to middle widtli 
is determined based on tlie corresponding immunity tran- 
sition points determined in step 2215 of all vectors. The 
approaches described above with reference to Figures 18, 
19 and 20 may be used to determine the curve for each 
vector. 

[0199] In step 2225, the accurate immunity transition point cor- 
responding to the middle width may be determined based 
on the NP curve determined in step 2220 for all vectors. 
The approaches described above with respect to Figures 
13B and 11 may be used to determine the accurate immu- 
nity transition point similar to step 2215. 

[0200] In step 2230, the (approximate) NIC is determined for all 
the vectors based on the three immunity transition points 
corresponding to narrowest, middle and widest widths us- 
ing techniques such as curve fitting as noted above. In 
step 2235, the vector corresponding to the lowest noise 
immunity curve is set as the worst case vector. Due to the 
use of the data corresponding to only the worst case vec- 
tor, computations are avoided for other vectors to at least 
to some degree. The remaining steps of Figure 22 are 
performed only on worst case vector. 

[0201] In step 2240, the start point corresponding to each of re- 



maining widths is determined according to tlie corre- 
sponding lieiglit on NIC. Due to tlie use of the data from 
NIC (for the worst case vector in the described embodi- 
ment), the start point may be quickly placed in the transi- 
tion region, thereby reducing the computational require- 
ments. 

[0202] In step 2250, the NP curve corresponding to each of the 

remaining width is determined using the start point deter- 
mined in step 2240. The determination may be performed 
using the approaches of Figure 18. Due to setting of out- 
put peaks to Vdd or 0 corresponding to all the input 
glitches that cause Vdd or 0 in the output peak, computa- 
tional requirements are reduced. 

[0203] In step 2260, the immunity transition points correspond- 
ing to each of the remaining widths are determined using 
the NP curve of the corresponding width of the worst case 
vector. The determination may be performed using the 
approaches of Figure 15. Due to the narrow search range 
obtained as described above with reference to Figure 15, 
computational requirements are also reduced accordingly. 

[0204] In step 2290, an integrated circuit is analyzed using the 
NP curves and NICs. The flowchart ends in step 2299. 
From the above, it may be appreciated that the computa- 



tional requirements are reduced by several techniques. 
Tlie description is continued with respect to reducing 
computational requirements in the characterization of se- 
quential elements. 

[0205] IX, Sequential Elements 

[0206] IX.A. General Introduction 

[0207] Noise immunity analysis often needs to be performed for 
various heights and widths of input glitches applied to the 
input pins of sequential cells. A noise immunity curve 
(NIC) for sequential cells needs to indicate the minimum 
height of an input glitch (of a specific width) which would 
cause a failure, i.e., generate an output glitch of a pre- 
specified height (e.g., 10% of Vdd) for any timing of a 
clock signal (that would determine the time instance at 
which the sequential element would change state). Such 
an output glitch would be generated when the setup and 
hold times are met with the relative timing of the input 
glitch and the clock signal. Thus, noise immunity analysis 
is performed on pins that are constrained by a clock sig- 
nal, for example, data pin. 

[0208] Thus, for input glitch 2310 of Figure 23 with a specific 

height and width, the clock transition point 2301 is varied 



to various points of interest to generate corresponding 
clocl< signals, sliown as dotted lines. The transition point 
is varied in sweep region 2320, which extends by a nega- 
tive setup time ahead of the expected start time of input 
glitch 2310 and by a negative hold time after the ex- 
pected end of the input glitch 2310. Generating such 
multiple clock signals is referred to as 'clock sweeping'. 

[0209] In one prior embodiment, the clock sweeping is per- 
formed for each combination of height and width of inter- 
est. Such an approach may require a large number of sim- 
ulations, and thus computation resources. For example, 
for an input noise glitch of width lOOOPS of height as- 
suming negative setup and negative hold margin of 20% 
each, total width for which clock sweeping will be required 
= 1400PS. If clock is swept at the step of lOPS, then total 
number of spice simulations equals 140. Thus 140 spice 
simulations may be performed for each combination of 
height and width of interest. The number of combinations 
of interest may be reduced by employing binary search 
type techniques within the 0 to 2 * Vdd range. 

[0210] Various aspects of the present invention reduce the num- 
ber of simulations further as described below. 

[021 1 ] IX.B. Avoiding Sweeps 



[0212] Figure 24 is a flowchart illustrating the manner in which 
sweeps can be avoided for some heights of input glitches 
according to an aspect of the present invention. The 
flowchart is described with reference to Figure 25. The 
flowchart starts in step 2401 in which control immediately 
passes to step 2410. 

[0213] In step 2410, a first value indicating the height of an input 
glitch is received. In step 2430, a simulation is performed 
on the (design of) cell with a clock signal providing large 
setup and hold times with respect to an input transition 
waveform having height equaling the first value. The large 
setup and hold times are often referred to as infinite 
setup and hold times. In general, the length/ duration of 
setup and hold times are chosen to ensure that the prop- 
agation of input glitch to the output is not limited by 
setup and hold time considerations as described below 
with reference to Figure 25. 

[0214] Figure 25 is a timing diagram illustrating the manner in 
which step 2430 may be performed. As shown there, in- 
put transition waveform 2510 is held at the height of the 
input glitch (Gh) for a long duration such that both setup 
time 2511 and hold time 2512 are substantially high/ 
infinite (with reference to clock signal 2520). The height 



of input transition waveform 2510 equals tlie lieiglit of in- 
put glitcli 2530, sliown as dotted line to the extent there 
is no overlap with input transition waveform 2510. Due to 
the high values of setup time and hold time, the latching 
of high level of input transition waveform 2510 is not lim- 
ited by setup time and hold time considerations. 

[0215] Continuing with reference to Figure 24, in step 2450, a 

determination is made as to whether a transition occurs at 
the output for the simulation of the input glitch. Control 
passes to step 2460 if a transition occurs and to step 
2480 otherwise. 

[0216] In step 2460, the clock signal is swept in the desired 

sweep region to determine whether a failure result occurs 
in at least one of the simulations. As with noise immunity 
analysis described above for combinatorial cells, a failure 
result is deemed to be obtained if the height of the output 
glitch exceeds 10%of Vdd (supply voltage/maximum 
height of the output glitch). Otherwise, a success result is 
deemed to be obtained. The sweeping may be stopped 
upon the occurrence of first failure result. Control then 
passes to step 2499, including in situations when no fail- 
ure results are obtained. 

[0217] In step 2480, it is concluded that input glitches having 



height of the first value (are not harmful) would not effect 
the cell. It may be inferred that the input glitches having 
the height of the first value with any width are also not 
harmful and can be ignored for performing simulation. 
Control passes to step 2499, in which the method ends. 

[0218] Thus, by using an extra simulation, sweeping (140 simu- 
lations in the above example) may be avoided for each 
width of the input glitch having a height of the first value, 
thereby reducing the computational requirements. The 
description is continued with respect to the manner in 
which noise immunity characterization can be performed 
for sequential elements according to various aspects of 
the present invention. 

[021 9] jx,c. Noise Immunity Characterization for Sequential Elements 

[0220] Figure 26 is a flowchart illustrating the manner in which 
noise immunity curves can be generated (for a specific 
width of input glitch) for sequential elements according to 
an aspect of the present invention. The flowchart starts in 
step 2601 in which control immediately passes to step 
2610. 

[0221] In step 2610, a search range is formed for a height of an 
input glitch. The search range may be formed based on 
the techniques described above with reference to Figures 



11 and 14A. In step 2620, a point in tlie search range for 
a present iteration is determined. Assuming binary searcli 
tecliniques are used, the middle point between an upper 
bound and lower bound of the search range is set as the 
point. 

[0222] In step 2630, a determination is made as to whether there 
is a failure result for the point. The approach of Figures 
24 and 25 can be used to make such a determina- 
tion. described, such a determination can be potentially 
made with one simulation. 

[0223] In step 2640 a determination is made whether to continue 
the search. In general, once the point with a failure result 
is determined with a desired resolution/accuracy, the 
search can be terminated. Control passes to step 2650 if 
the search is to be continued and to step 2699 otherwise. 
In step 2650, the next point for simulation is determined 
according to binary search approach and control passes to 
step 2630. The method ends in step 2699. 

[0224] By using narrow search ranges in step 2610 (in addition to 
the techniques of Figures 24 and 25), the computational 
requirements in generating NIC may be reduced. The de- 
scription is continued with respect to embodiments of 
digital processing system implemented substantially in 



the form of software. 

[0225] X. Software Implementation 

[0226] Figure 27 is a block diagram illustrating the details of dig- 
ital processing system 2700 implemented substantially in 
the form of software in an embodiment of the present in- 
vention. System 2700 may contain one or more proces- 
sors such as central processing unit (CPU) 2710, random 
access memory (RAM) 2720, secondary memory 2730, 
graphics controller 2760, display unit 2770, network in- 
terface 2780, and input interface 2790. All the compo- 
nents except display unit 2770 may communicate with 
each other over communication path 2750, which may 
contain several buses as is well known in the relevant arts. 
The components of Figure 27 are described below in fur- 
ther detail. 

[0227] CPU 2710 may execute instructions stored in RAM 2720 
to provide several features of the present invention. CPU 
2710 may contain multiple processing units, with each 
processing unit potentially being designed for a specific 
task. Alternatively, CPU 2710 may contain only a single 
processing unit. RAM 2720 may receive instructions from 
secondary memory 2730 using communication path 2750. 

[0228] Graphics controller 2760 generates display signals (e.g.. 



in RGB format) to display unit 570 based on data/ 
instructions received from CPU 2710. Display unit 2770 
contains a display screen to display the images defined by 
the display signals. Input interface 2790 may correspond 
to a key.board and/or mouse, and generally enables a 
user to provide inputs. Network interface 2780 enables 
some of the inputs (and outputs) to be provided on a net- 
work. In general, display unit 2770, input interface 2790 
and network interface 2780 enable a user to check the 
status of various characterization tasks, and to perform 
chip-level analysis of integrated circuits. 
[0229] Secondary memory 2730 may contain hard drive 2735, 
flash memory 2736 and removable storage drive 2737. 
Secondary memory 2730 may store the data and software 
instructions (e.g., to characterize a parameter), which en- 
able system 2700 to provide several features in accor- 
dance with the present invention. Some or all of the data 
and instructions may be provided on removable storage 
unit 2740, and the data and instructions may be read and 
provided by removable storage drive 2737 to CPU 2710. 
Floppy drive, magnetic tape drive, CD_ROM drive, DVD 
Drive, Flash memory, removable memory chip (PCMCIA 
Card, EPROM) are examples of such removable storage 



drive 2737. 

[0230] Removable storage unit 2740 may be implemented using 
medium and storage formatcompatible with removable 
storage drive 2737 such that removable storage drive 
2737 can read the data and instructions. Thus, removable 
storage unit 2740 includes a computer readable storage 
medium having stored therein computer software and/or 
data. 

[0231] In this document, the term "computer program product"is 
used to generally refer toremovable storage unit 2740 or 
hard disk installed in hard drive 2735. These computer 
program products are means for providing software to 
system 2700. CPU 2710 may retrieve the software in- 
structions, and execute the instructions to provide various 
features of the present invention as described above. 

[0232] XI. Conclusion 

[0233] While various embodiments of the present invention have 
been described above, it should be understood that they 
have been presented byway of example only, and not 
limitation. Thus, the breadth and scope of the present in- 
vention should not be limited by any of the above de- 
scribed exemplary embodiments, but should be defined 
only in accordance with the following claims and their 



equivalents. 



